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TITLE OF THE INVENTION 

METHOD AND SYSTEM FOR SENSING AND ANALYZING A WAVEFRONT 

OF AN OPTICALLY TRANSMISSIVE SYSTEM 

BACKGROUND AND SUMMARY OF THE INVENTION 

[00011 Technical Field. 

[0002] This invention pertains to the field of wavefront sensing methods and devices, and 
more particularly, wavefront sensing methods and devices used to measure the optical quality 
of an optically transmissive system or device, for example, an optical component such as a 
lens. 

[0003] Description. 

[0004] A Ught wavefront may be defined as the virtual surface delimited by all possible 
rays having an equal optical path length from a spatially coherent source. For example, the 
wavefront of light emanating from a point light source is a sphere (or a partial sphere where 
light from the point light source is lunited to emission along a small range of angles). 
Meanwhile, the wavefront created by a collimatiag lens mounted at a location one focal 
length away from a point source is a plane. A wavefront may be planar, spherical, or have 
some arbitrary shape dictated by othCT elements of an optical system through which the ligjit 
is transmitted or reflected. 

[00051 A wavefront analysis system, including a wavefront sensor, may be used to measure 
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characteristics of an optically transmissive system by detecting the wavefront of light 
emerging from the system and comparing it to some expected ideal wavefront (planar, 
spherical, etc.). The optically transmissive system might be a single component or may be 
very complex, such as a transmissive optics system, (e.g., a coUimating lens; ophthalmic 
lens). The differences between the ideal expected wavefront and the actual measured 
wavefront are caused by optical aberrations of the system under test (SUT). 
[0006] A number of different wavefront sensors and associated methods are known. 
Among these are interferometers and the Shack-Hartmann wavefront sensor. Each of these 
systems will be described briefly below. A more detailed discussion of wavefront sensing 
techniques may be found in "Introduction to Wavefront Sensors," 1995, Joseph M. Geary, 
SPIE Press. 

[0007] Interferometers. 

[0008] An interferometer is an instrument that uses interference of light waves to detect the 
relative wavefront difference between a test light beam and a reference beam. Interferometric 
methods of sensing a wavefront are highly sensitive but very limited in dynamic range. A 
typical interferometer can only directly measure optical path differences of less than one 
wavelength - a 27C phase ambiguity exists beyond the one wavelength point. If the optical 
path difference is greater than one wavelength, then the correct phase difference is often 
inferred computationally using phase imwrapping techniques. However, real optical 
configurations can be constmcted where these techniques are likely to fail. Other limitations 

« 

of interferometric techniques include the necessity of relative stability of the reference and 
test beam paths. This means that any vibration in the test instrument leads to a degradation of 
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the measurement accxiracy. 

[00091 Shack-Hartmann Wavefront Sensors. 

[00010] A Shack-Hartmann wavefront sensor is a device that uses the fact that light travels 
in a straight Une, to measure the wavefiront of Ught. FIG, 2 shows a basic configuration of a 
Shack-Hartmann wavefront sensor 200. The Shack-Hartmann wavefront sensor 180 
comprises a lenslet array 1 82 that breaks an incoming beam into multiple focal spots 1 88 
falhng on an optical detector 184. Typically, the optical detector 184 comprises a pixel array, 
for example, a charge-coupled device (CCD) camera. By sensing the positions of the focal 
spots 188, the propagation vector of the sampled Ught can be calculated for each lenslet of the 
lenslet array 182. The wavefront can be reconstructed from these vectors. 
[00011] However, Shack-Hartmann wavefront sensors have a finite dynamic range 
determined by the need to associate a specific focal spot to the lenslet it represents. A typical 
methodology for accomplishing this is to divide the detector surface into regions (called 
"Areas Of Interest" [AOIs]) where the focal spot for a given lenslet is expected to fall. If the 
wavefront is sufficiently aberrated to cause the focal spot to fall outside this region, or not be 
formed at all, the wavefront is said to be out of the dynamic range of the sensor. FIG. 3 
shows an example of a Shack-Hartmann wavefront sensor 300 in an out-of-range condition. 
[00012] In practice Shack-Hartmann wavefront sensors have a much greater dynamic range 
than interferometric sensors. This range may be tens to hundreds of waves of optical path 
difference. However, this dynamic range is still insufficient to characterize many real optics. 
[00013] Other Wavefront Sensing Technologies. 

[00014] Other sensors such as the Moire Deflectometer have a higher dynamic range, but 
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lack the sensitivity necessary for accurate measurement of most transmissive optical 
elements. 

[00015] Both optical and computational methods have been used to extend the dynamic 

range of wavefront sensing devices. Some example computational methods include Spot 

« 

Tracking, Phase Unwrapping, and Angular Spectrum Propagator Reconstruction. 
[000161 Spot Tracking. 

[00017] This method extends the dynamic range of the Shack-Hartmaim wavefront sensor in 
the case where the wavefront being measured starts out within range and then drifts out of 
range over a period of time. This case exists for many optical configurations where a lens 
moves within the optical setup or a component changes optical characteristics due to some 
cause such as material heating or deformation. Spot tracking is accomplished by comparing 
current positions of focal spots to positions recorded in a previous frame. The previous 
positions are used as a starting point for locating the spots after an incremental movement. 
As long as the frames are taken frequently enough, then it is computationally simple to keep 
track of them. This technique has been known since at least 1993 (A. Wirth, A Jankovics, F. 
Landers, C. Baird, and T. Berkopec, 'Tinal report on the testing of the CIRS telescopes using 
the Hartmann technique," Tech. Rep. NAS-31786, Task 013 (Adaptive Optics Associates, 
Cambridge, MA 1993)). A limitation to this approach is that the incident wavefront must 
start out within range. 

[00018] Phase Unwrapping. ; 

[00019] In this technique the focal spot to lenslet mapping is inferred using techniques 
similar to those used in interferometry. This technique is described in "Dynamic range 
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expansion of a Shack-Hartmann sensor by use of a modified ujiwrapping algorithm," by J. 

Pfimd, N. Lindlein, and J. Schwider, Optical Society of America, 1998. 

> 

[00020] Angular Spectrum Propagator Reconstruction. 

[00021] Described by "Algorithm to increase the largest aberration that can be reconstructed 

« 

from Hartmann sensor measurements," by M. Roggemann and T. Shulz, Applied Optics, Vol 
37, No 20, 1998, this technique is computationally expensive and therefore inappropriate for 
many measurement applications. 

[00022] While there are means for tracking and adjusting the positions of these AOIs (as 
described previoxisly), the simplest, most robust calculations are achieved for the case where a 
single mapping of lenslets onto the pixels can be maintained. For example, U.S. patent no. 
5,825,476 discloses a method that uses a naissing focal spot to identify the central AOI, and 
then tracks all the other focal spots using this missing data. However, if there is a speck of 
dust on the part under test, this easily fools the identification of this missing spot, leading to 
inaccurate results. 

[00023] U.S. Patent No. 6,550,917, the entire contents of wluch are hereby incorporated by 
reference for all purposes as if fiilly set forth herein, discloses a means for extending the 
dynamic range of a sensor by adjusting the spherical radius of curvature of a reference sphere 
to match the effective defocus of the optical system imder test (in that case, an eye). While a 
similar scheme can be applied to testing transmissive optics, the requirements for testing an 
eye en vivo are significantly different from those of measuring a fixed lens or optical element. 

For an intraocular lens, the focal length in air can be as short as 10 mm, necessitating the use 

of a different optical testing method. 
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[00024] Without some tracking scheme, the focal spot moves either completely or partially 
out of the assigned AOL Thus the pixels do not accurately represent the focal spot position 
information. If inaccurate data is used to set the position of an adaptive focal element, then 
the system may not converge to the correct position, leading to an inaccurate measurement or 
slow convergence. In a typical embodiment, there may be thousands of potential focal spots. 
Thus identifying poor focal spots, or focal spots that have wandered outside of their correct 
AOIs, is useful for maintaining the accuracy of the subsequent calculations. 
[00025] Accordingly, it would be desirable to provide a system and method for extending the 
dynamic range of wavefront sensing devices in transmissive optics metrology. It would also 
be desirable to provide a system and method for sensing and analyzing the wavefiont of light 
passing through an optically transmissive system with enhanced dynamic range. It would be 
further desirable to provide such a system and method which overcomes one or more 
disadvantages of the prior art. 

[00026] The present invention comprises a system and method for sensing and analyzmg a 
wavefrpnt of an optically transmissive system. By analyzing the wavefront, the system and 
method may ascertain desired parameters of the optically tranismissive system. Li particular, 
for example, when the optically transmissive system is a lens. such as a contact lens, the 
system and method may accurately determine the focal length of the lens by sensing and 

■ 

analyzing the wavefront of light passing through the lens. Furthermore, once the focal length 
has been determined, the system and method may eliminate the lower order focal length term 

from the analysis, and determine higher order aberrations in the lens. 

[00027] In one aspect of the invention, a system for sensing a wavefront of light passed 



Docket No. WFS.OllC 
-PATENT-- 

through an optical device, comprises: a device under test (DUT) holder adapted to hold the 
optical device; a point light source adapted to provide light to the optical device; a movable 
platfomi adapted to move the point light source along an optical axis of the system; a first 
lens adapted to receive and pass therethrough a light beam from the optical device; a range- 
limiting aperture adapted to receive and pass therethrough at least a portion of the light beam 
from the first lens; a second lens adapted to receive and pass therethrough the portion of the 
light beam from the range-limiting aperture; a Shack-Hartmann wavefront sensor adapted to 
receive the portion of the light beam from tihe second lens and to produce therefrom 
wavefront data; and a processor adapted to receive the wavefront data from the wavefront 
sensor and to control movement of the movable platform to move the point light source to a 
location about one focal length away from the optical device. 

[00028] In another aspect of the invention, a method of measuring a wavefront of light from 
an optically transmissive device, comprises: locating a light source a first distance from the 
optically transmissive device; passing light from the light source through the optically 
transmissive device; imaging at least a portion of the light pa$sed through the optically 
transmissive device; sensing a wavefront of the imaged light to produce therefrom wavefront 
data; and adjusting a location of the light source witii respect to the optically transmissive 
device to substantially maximize a degree of collimation of the light passed through the 
optically transmissive device. 

[00029] In yet another aspect of the invention, a system for measuring a wavefront passed 
through an optically transmissive device, comprises: a light source disposed on a first side of 
an optically transmissive device; a wavefront sensor disposed on a second side of an optically 
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transmissive device; a relay imaging system disposed between the optically transmissive 
device and the wavefront sensor; and means for adjusting a distance between the light source 
and the optically transmissive device. . ■ 

[00030] In still another aspect of the invention, a method of determining a focal length of a 
lens comprises: (1) locating a light source a first distance from the optically transmissive 

4 

device; (2) passing light from the light source through the opticklly transmissive device; (3) 
imaging at least a portion of the ligjit passed through the optically transmissive device; (4) 
sensing a wavefront of the imaged light; adjusting a location of the light source with respect 
to the optically transmissive device to substantially maximize a degree of collimation of the 
light passed througji the optically transmissive device; (5) moving the light source by a 
distance Xi from the location that substantially maximizes the degree of collimation of the 
light passed through the optically transmissive device, where i=(l, N); (6) sensing a 
wavefront of the imaged light; (7) calculating a radius of curvature value of the sensed 
wavefront; (8) repeating the steps (5) through (7) N-1 times (where N is an integer) where the 
value of Xi is changed each time the steps (5) through (7) are repeated; and (9) calculating the 
focal length of the lens from the N radii of curvature values calculated in the steps (5) though 
(8). i 

BIOEF DESCRIPTION OF THE DRAWINGS 
[00031] FIG. 1 is a functional diagram of one embodiment of a wavefront analysis system 

■ 

according to one or more aspects of the invention; 
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[00032] FIG. 2 shows a basic arrangement of a Shack-Hartmaim wavefiront sensor; 
[00033] FIG. 3 shows an out-of-range condition for a Shack-Hartmann wavefiront sensor; 
[00034] FIG. 4 shows an example area of interest (AOI) for a case where the focal spot is 
well within range of a lenslet; 

[0003S] FIG. 5 shows an AOI that has light from two different lenslets; 

[00036] FIG. 6 shows a case where there is sufficient tilt of a light beam from a system under 

test that a range-limiting aperture partially occludes the focal spot; 

[00037] FIG. 7 shows the pattern of focal spots for a system such as that shown in FIG. 1 

when a light beam emerging from the system under tests is highly collimated; 

[00038] FIG. 8 shows the pattem of focal spots for a system such as that shown in FIG. 1 

when a light beam emerging from the system under tests is less coUimated than in FIG. 8; 

[00039] FIGs. 9A-C illustrate operation of a range-limiting aperture in conjunction with a 

wavefront sensor imder three different conditions; and 

[00040] FIG. 10 is a diagram for illustrating a method of determining a focal length of a lens. 

DETAILED DESCRIPTION 

[00041] FIG. 1 shows a functional diagram of an embodiment of wavefront analysis system 
100. The wavefront analysis system 100 includes a point light source 1 10, a movable platform 
or stage 120, a system under test (SUT) holder 130, an objective lens 150, an aperture 160, an 
eyepiece lens 170, a wavefront sensor 180, and a processor 190. 

[00042] An exemplary device imder test (SUT) 140 is also shown in FIG. 1, comprising an 
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ophthalmic lens. The SUT 140 may be any optically transmissiye system or device, including 
but not limited to an optical component. Specific examples include a collimating lens, a 
contact lens, an intra-ocular lens, plate glass, a windshield, etc., or a collection of optical 
devices. 

[00043] The point Ught source 1 10 may be a laser or laser diode, but, for example, it may 
also be an incandescent white light source coupled through a fiber or pin hole to form a point 
source, or any similar arrangement. Such a configuration has an advantage that, by adding 
color filters, chromatic aberrations of the SUT 140 can be measured. 

[00044] In some cases it may be desirable to replace the point Ught source with another light 
source of known divergence. This would allow the measurenient to be scaled to different 
ranges of radius or curvature. To this end additional optical elements would be added in front 
of the point source to provide a means for setting the divergence to a known value. 
[00045] The movable platform 120 is a mechanism that allows the location of the point light 
source 1 10 to be moved relative to the SUT holder 130. The movable platform 120 may be 
an electro-mechanical stage, a stage with a micrometer, or a similar contrivance. The 
movable platform 120 may be manually adjusted. However, beneficially the movable 
platform 120 is adapted to be controlled by a feedback signal to automatically move the point 
light source 110 to a desired location with respect to the SUT 140, as wiU be described in 
more detail below. For example, the movable platform 120 may include a motor that 
transports it using rollers, wheels, along a rail, etc. Also beneficially, the movable platform 
120 maybe attached to a position digitizer that produces location data indicating a relative 
distance between the movable platform 120 and the SUT 140. 
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[00046] As the movable platform 120 with the point light source 1 1 0 is moved, it may be 

• « ■ 

desirable to adjust the amount of light either eroitted from the point light source 110, or 
detected by the wavejBront sensor 180 (through integration time or gain control) so that the 
wavefront sensor 180 operates at an appropriate Ught level. This adjustment can be 

a 

performed automatically in a closed loop control circuit, or manually if only small or 
occasional adjustments are needed. 

4 

[000471 The objective lens 150 and an eyepiece lens 170 form a Keplerian telescope and 
together form an imaging system, i.e., a relay imaging system. 

[00048] Beneficially, the wavefront sensor 180 is a Shack-Hartmann wavefront sensor as 
shown in FIG. 1, including a lenslet array 182 producing a corresponding plurality of light 
spots 188 firom the received light beam, and an optical detector comprising a pixel array (e.g., 
a CCD camera) 184 upon which the light spots 188 are focused. For each lenslet in the 
lenslet array 182, a subset of the pixel array 184 defmes an "Area-of-Interest" (AOI). For 
example, where there are 1200 lenslets in the lenslet array 182, and where the pixel array 184 
comprises 640 by 480 pixels, then the AOI assigned for each lenslet may be a square region 
of 16 by 16 pixels (256 pixels). Accordingly, characteristics of the wavefiront of ttie received 
Ught beam may be ascertained by examining where the light spot 188 is projected in each 
AOI of the pixel array. 

[00049] The aperture 160 is located at the inter-lens focal point of the imaging system for the 
selected focal power of the lenses. The aperture 160 is arranged such that it has appropriate 
size and shape to match the size and shape of the AOI after imaging througji the optical 
system consisting of lens 170 and lenslet array 182 of the wavefront sensor 180. The shape of 



11 



Docket No. WFS.OllC 
--PATENT-- 

aperture 160 maybe chosen as square, hexagonal, round or other convenient shape so long as 
it images to a region substantially the sanie size as the AOI for each lenslet in the lenslet array 
182 of the wavefiront sensor 180. Accordingly, the aperture 160 will only transmit light for a 
small range of input radii-of-curvature ligjit beams. Typically, the size of the aperture 160 is 
such so that the system 100 is transmissive only for light beaitis that are very nearly 
collimated. 

[00050] The processor 190 is adapted to receive and processes wavefiront sensing data from 
the wavefront sensor 180. Beneficially, the processor 190 may process the wavefront sensing 
data and generate therefrom a feedback signal which it provides to the movable platform 120, 
as described in more detail below. For example, the feedback signal may be proportional to 
the curvature of the wavefront of the light beam received by the wavefront sensor 180. The 
feedback signal may be produced by an algorithm, executed by the processor 190, that 
determines the relative curvature of the wavefront by comparing the separation of Ught spots 
on the pixel array of the wavefront sensor 180. 

[00051] The operation of the wavefront analysis system 100 will now be described. 
[00052] Light from the point light source 1 10 is received by and passes through the SUT 
140, producing a light beam. The light beam from the SUT 140 is received by the imaging 
system comprising the objective lens 150 and the eyepiece len§ 170. The imaging system 
sizes the light beam appropriately for the wavefront sensor 180 and performs relay imaging, 
allowing the wavefront measurement to be made at the appropriate plane. The imaging 
system also provides an internal focus point for the aperture 160. 

[00053] At least a portion of the light beam passes through the aperture 160 to the eyepiece 
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lens 170 and therefrom to the wavefront sensor 180. 

[00054] The lenslet array 1 82 and second lens 170 of the relay imaging telescope form an 
imaging system between detector array 184 and the aperture 160. Because the lens 170 is a 
single lens, with an aperture much greater than the individual, lenslets of the lenslet array 182, 
a common point at the middle of aperture 1 60 will form an image on the detector array 184 
directly behind each individual lenslet. That is, the image wo^ld form on the center of each 
AOL Variations in the local wavefront slope at varioxis locations of the SUT 140 result in a 
shift of the position of the focal spot at the aperture 160, and, since the system 170 and 182 
act as an imaging system, also result in a position shift of the focal spots at the detector array 
184. It is exactly this shift that is detected by the wavefront sensor 180 as a measure of the 
local wavefront slope. 

[000551 The aperture 160 ensures that the wavefront received by the wavefront sensor 180 
remains within the initial dynamic range of the wavefront sensor 180. For the case where 

local wavefront slope on the SUT 140 is irregular (either because of optical errors or due to 

t 

fact that the source 110 and device imder test 130 are not arranged to produce a coUimated 
beam, that is, a planar wavefront), then the focal spots 188 will move across the detector array 
184 of the wavefront sensor 1 80 in a pattern that corresponds to the irregular wavefront 
slopes. If the local wavefront slope is too great (the out-of-range condition), it will be clipped 
by the aperture 160 and hence not be incident upon the detector array 184. Thus, in this 
arrangement, the aperture 160 removes any light that otherwise would cause an out-of-range 
condition for the lenslet array 182 of the wavefront sensor 180. 

[00056] Thus the aperture 160 serves as a range-limiting aperture, thereby limiting the 
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dynamic range of the incident light so that the wavefiront sensor 180 never receives light that 
would otherwise cause an out-of-range condition. Typically, the hght incident at the aperture 
160 creates a focal spot that is less than ten (10) times the diameter of the aperture 160 itself, 
although this depends upon the particular system design. 

[00057] When the point Ught source is located far away firom the optimal position, the Ught 
passing througji the objective lens 150 will include components that are severely off-axis and 
which arrive at the plane the range-limiting aperture 160 spaced fairly far away from the 
internal focal point of the imaging system. Thus, these components are blocked by the 
range-limiting aperture 160. This insures that any out-of-range condition of the wavefront 
sensor 180 catmot occur. 

[000581 FlGs. 9A-C illustrate operation of the range-limiting aperture 160 and the wavefront 

■ 

sensor 180 under three different conditions. In FIG. 9A, all of the light enters through the 
center of the range-limiting aperture 160. FIG. 9B illustrates the case where the light enters 
near the edge of the range-limiting aperture 160. In that case, the light spots 188 all still 
resolve within a corresponding AOI at the detector array 184. FIG. 9C illustrates the case 
where the light enters at the edge of the range-limiting aperture 160. 
[00059] Advantageously, the point light source 1 10 is configured so that the distance 
between the point light source 110 and the SUT 140 may be adjusted. Beneficially, the point 
light source is attached to the movable platform 120 and thereby may be moved along the 
optical axis of the wavefront analysis system 100. 

[00060] Beneficially, when the SUT 140 is a contact lens or intra-ocular lens, for example, 
the wavefront analysis system 100 may be ernployed to determine the focal length of the SUT 
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140. 

[00061] In the wavefiront analysis system 100, the light beam, provided to the wavefiront 
sensor 180 is more or less collimated depending upon the rela^ye location of the point light 
source 1 10 with respect to the SUT 140. When the point light source 1 10 is located one focal 
length away from the SUT 140, then the light beam from the SUT 140 will be optimally 
collimated. During the initial set-up, the wavefront sensor 180 and processor 190 determine 
the degree of collimation of the light beam produced by the point light soxirce 110 and the 
SUT 140. Then, the point light source 1 10 is moved by the movable platform 120 until it is 
located at a position with respect to the SUT 140 that produces an optimal or near-optimal 
degree of collimation of the light beam emerging from the SUT 140. 
[00062] Although the movable platform 120 may be moved manually, beneficially an 
electro-mechanical movable stage 120 is employed which operates under the control of the 
processor 190. In that case, the processor 190 may operate with the movable stage 120 in a 
feedback loop to move the point Ught source 1 10 to a location along the optical axis of the 
system 100 where the ligjit beam from the SUT 140 is optimally or nearly optimally 
collimated. And the position of the point light source 110 may be continuously adjusted as 
necessary to accoimt for any vibration or changes in temperature that may change the optimal 
location for the point light source 1 10. 

[00063] The distance between the point light source 1 10 and the SUT 140 that produces 
optimal collimation of the light beam may be interpreted as the focal point (focal term) of the 
SUT 140. 

[00064] In theory, then, once the point light source 1 10 is located where it produces optimal 
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or near-optimal coUimation of the Eght beam from the SUT 140, in theory the focal length of 
the SUT 140 could be determined my measxiring the distance between the point light source 
110 and the SUT 140. This may be done in several ways, either jnanually or in an automated 
fashion. Beneficially, the location of the point light source 110 with respect to the SUT 140 
may be automatically detected from a position digitizer 125 by the processor 190. 
[00065] However, such a direct measurement approach for determining the focal length of 
the SUT 140 suffers from inaccuracies due to a number of factors, including, e.g., the finite 
thickness of the contact lens/intra-ocular lens itself. In many applications, tiie inaccuracy of 
tiie focal length value produced by the direct measurement approach is too great. 
[00066] Accordingly, another method has been provided to employ the wavefiront analysis 
system 100 to determine the focal length of the SUT 140. This method will hereaft^ be 
described with respect to FIG. 10, which shows a lens 1000, which may be the SUT 140 in 
FIG. 1. In FIG. 10: "fi" indicates the location where a point light source should be placed 
such that the lens 1000 will output a fully coUimated light (i.e., fi is located one focal length 
away from lens 1000); "x" indicates a distance between the location "fi" and the location 
where a point light source is located; "f2** indicates a location one focal length away from 
lens 1000 on the other side of the lens 1000 from the point light source; "e" indicates a 

■ 

distance between the location "fi" and the location where a wavefront sensor (such as the 
wavefront sensor 180 of FIG. 1) is located; *¥' indicates the radius of curvature of a ligjit 
beam passing through lens 1 000 when illuminated by a point ligjit source placed at the 
location "x'*; and 'y ' indicates a distance between the location "fa" and the location where 
the light beam passing through the lens 1000 converges when illuminated by a point light 
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source placed at the location "x." 

t 

[00067] In FIG. 10, when a point light soiirce is located at a distance "x" from the first focal 
point "fi " and a wavefront sensor is located at a distance "e" from ttie second focal point 
"fi," then according to Newton's thin fihn formula as shown in Equation (1) below. 

[000681 1) x(rfe) = -f^ 

[00069] where "f * = the focal length of the lens 1000, and, as noted above, *V = the radius 
of curvature of the Ught beam passing through the lens 1000 as measured by the wavefront 
sensor. 

[00070] So, if 'V and "e" are known, one can solve Equation (1) for "f by measuring "r," 
the radius of curvature of the light beam passing through the lens 1000, with a wavefront 
sensor. Furthermore, even if "e" is not known, one can solve for the two unknown variables 
and "e'* by solving Equation (1) for two different values of "x" (two different point light 
source locations) and "r.** 

[00071] Accordingly, a method for determining the focal length "f of the lens 1000 with the 
wavefront analysis system 100 of FIG. 1 proceeds as follows. , 

[00072] First, the wavefront sensor 1 80 is placed an unknown distance ("e'*) downstream 
from the SUT 140 (e.g., the lens 1 000). 

[00073] Next, the location of "fi" is determined by moying the point light source 1 10 via the 
movable platform 120 until it is located at a position with respect to the SUT 140 that 
produces an optimal or near-optimal degree of coUimatipn of the light beam emerging from 
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the SUT 140, as detennined by the wavefro^t sensor 180. Specifically, data measvired by the 

wavefront setisor 180 is used to calculate a radius of curvature *¥' of the light beam, and the 
movable platform 120 is adjusted until *V* is maximized (goes toward an infinite radius of 
curvature, which would indicate a perfectly coUimated light beam). The radius of curvature 
"r" of the light beam may be determined firom the data detected by the wavefront sensor 180 
according to a variety of algorithms. For example, the wavefront data may be used to 
compute the coefficients of the so-called Zemike polynomials. In that case, the radius of 
curvature of the light beam may be determined from equation (2) below. 



[000741 2) "r" = / (4* Z2.0) 



[00075] where R = radius of the lens 1000, and Zi.o is the coefficient of the (2,0) term of the 
Zemike polynomial. 

[000761 Furthermore, an exemplary algorithm for calculating the radius of curvature "r" 
ftom the senses wavefront data is described in W.H. Southwell, '*Wave-front estimation from 
wave-front slope measurements," J. Opt. Soc. Am., Vol. 70, No. 8, pp. 998-1006 (Aug. 1980) 
*the Southwell paper," the entirety of which is hereby incorporated herein in its entirety for 
all purposes as if fully set forth herein. 

[000771 Once the point of optimum collimation is reached, the point Ught source 1 1 0 is 
located at "fi" in FIG. 10. At this point, the location "fi" may be automatically detected from 
the position digitizer 125 by the processor 1 90. 

[000781 Now that the location "fi" has been determined, the movable stage 120 may be 
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employed to move the point light source 110 a distance **x" from "ff* and the distance "x" 
may be automatically detected from the position digitizer 125 by the processor 190. As the 
point light source 110 moves away ^piq. the focal point "fi " the wavefront sensor 180 will 
record smaller and smaller radii of curvature "r." 

[00079] Returning again to Equation (1), with the variable "x" being measured by the 
position digitizer 125, and the variable "r*' being determined by the wavefront sensor 180, 
there remain only two unknown variables in Equation (1): "e" and "f * (which is the focal 
length; the value we wish to deteraiine). Therefore, if the movable stage 120 is employed to 
move the point light source 1 10 to two differmt distances "x" from the focal point "fi," and if 
the two different corresponding radii of curvature are determined by the wavefront sensor 
180, then one can simultaneously solve the Equation (1) "twice" for the two unknown 
variables. 

[000801 Although in principle the focal length "f * can be detemiined by simultaneously 

« 

solving Equation (1) twice for two different values or "x" and **r," employing only two 
measurements and solving only two equations may still limit the accuracy with which the 
focal length "f is determined. 

[00081] Accordingly, beneficially, the movable stage is employed to place the point Ught 
source 100 at many CTSl*') different distances "x" from the focal point "fi," where i = (1, N). 
For each distance "x " where the point light source 1 10 is located, the wavefront sensor 180 
deteraiines the corresponding radius of curvature "r " of the Ught beam, where i = (1, N)- 
Then, an algorithm is applied to the data to determine the '*best fit" values for "e'' and "f." 
Beneficially, the algorithm may be a least-squares fit algorithm, although other techniques 
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may be employed. 

[00082] Details of one exemplary least-squares best fit algorithm may be found in, 
"Benjamin J. Pemick et al., "Least-squares technique for determining principle plane location 
and focal length/' AppUed Optics, Vol. 26, No. 15, pp. 2938-39 (Aug. 1987), the entirety of 
which is hereby incorporated herein in its entirety for all purposes as if fully set forth herein. 
[00083] However, employing such a least-squares best fit algorithm may yield a result that is 
more inaccurate than is desired if the data point at "x' -0 (at the focal point "fi") is employed 
in the algorithm. For example, Tables 1 and 2 below show some experimental results that 
were obtained with a lens having a focal length of 1.00 ± 0.05 m. 

Table 1 



x,mm r, mm 

60 15782.46325 

50 18891.13451 

40 23960.75208 

r' 

30 32353.13494 

20. 50591.44026 

10 104688.7173 

0 -1&39049.364 
-10 -103184.3176 
-20 -50736.61175 
-30 -33758.93702 
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-40 



-25711.53712 



-50 



-20803.82136 



-60 



-17466.85276 



Table 2 



focal length "f ' 



0.963 ± 0.009 m 



offset distance "e" 



0.8 ± 100 m 



[00084] It can be seen that the computed focal length "f ' is inconsistent with the expected 
value of 1.00 ± 0.05 m. And while the value for the offset distance "e" is reasonable, the 
uncertainty (±100 m) suggests that this value is nebulous. 

[00085] Accordingly, it has been discovered that a more accurate determination of the focal 
length "f is produced when any data point at or very near to "x"=0 is removed from the data 
set prior to employing the least-mean-square fit algorithm. When the point at "x**=0 is 
excluded from the data set of Table 1, then the values shown below in Table 3 are obtained, 
instead of those in Table 2 . 



[00086] As an alternative to eliminating the points neax x=0, it may be appropriate to 



Table 3 



focal length "f * 



1.002 ± 0.002 m 



offset distance "e" 



0.8 ± 8 m 
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formulate the fitting problem in terms of the dioptric power of the optical system as measured 
by the wavejfront sensor 1 80. That is, instead of fitting to curves of "r" vs. "x'\ it would be 
better to fit to curves of "l/r^* vs. "x". In this case the fit is not nearly so nonlinear, and a 
much more stable solution can be obtained. 

[00087] Finally, in some cases there may be uncertainly as to the precise location "fi" where 
the coUimation of the light beam is maximized. In that case, beneficially a least-squares best 
fit algorithm is employed to solve for THREE xmknown variables C^x," **e," and "f 0 instead 
of just two. It has been foimd empirically that this often yields more accurate results than 
when "fl" and **x** are determined by the position digitizer 125 as described above. This may 
be applied to either the direct fit for the radius of curvature *Y', or for dioptric power "1/r". 
[00088] Meanwhile, real lenses are not ideally thin. In general, the lens has a firont and a 
back principal plane. For example, the back principal plane of the lens is located somewhere 
inside the material of the lens. Often, the location of the back principle plane of the lens 1000 
is also a parameter of interest, which it is desired to measure with the wavefront analysis 
system 100. 

■» 

[00089] In the above-described algorithms, the parameter "e" is also calculated. Meanwhile, 
the distance **K'* firom the wavefi-ont sensor 180 to the back principal plane of the lens 1000 



equals "e" + "f^." 

[00090] Accordingly, in one embodiment, the wavefront analysis system 100 includes a 
second camera to precisely determine a distance "J" firom the wavefiront sensor 180 to the 
physical back surface of the lens 1000. The second camera looks at the lens 1000 in a 
direction perpendicular to the optical axis. Knowing 'TK." firom the algorithm described 
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above, and knowing "J" from the second camera, one can determine the location of the back 
principle plane of the lens 1000 with respect to its. back physical siirface by subtracting "J" 

from 'x;' 

[00091] The wavefront analysis system 100 may also be employed to characterize any 
surface irregularities and other imperfections in the SUT 140. In that case, beneficially the 
basic focal length term of the SUT 140 is eliminated from the aberrations of the wavefront 
that are measured directly by the wavefront sensor 180 by positioning the point light source 
110 with respect to the SUT 140 to produce optimal collimation of the light beam detected by 
the wavefront sensor 180. Accordingly, after having eliminating the dominant lower-order 
focal length term of the SUT 140, the data from the wavefront sensor 180 may then be 
employed to detect the higher order or residual optical aberrations of wavefront passing 
through the SUT 140. Such higher-order aberration data may be employed to assess the 
degree of surface irregularities and imperfections in the SUT 140. Further calculations such 
as the point spread ftinction (PSF), niodulation transfer ftmction (MTF), or other parameters 
maybe derived from the wavefront data. This is important for some optical systems since 
standards based on other measurement methods (e.g. MTF) maybe used for comparison to 
the wavefront data. 

[000921 A principle advantage of this system 100 and the method described above is that by 
removing the focal term from the analysis, the effective dynamic range of the wavefront 
analysis system 100 may be greatly increased. Since there is a well-known trade-off between 
dynamic range and sensitivity in typical wavefront sensor technologies, the wavefront 
analysis system 100 and associated method discussed above may incorporate a more sensitive 
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(and more accurate) wavefront sensor 1 80. Meanwhile, the range-limiting aperture 1 60 
insures that there is no location of the point light source 110 that produces an out-of-range 
condition for the waveftont sensor 180 at any point on the detector 184. This enables the 
processor 190 to adjust the distance between the point light source 110 and tiie SUT 140 (e,g., 
by moving the naovable platform 120), through the feedback loop in a stable manner that 
ensures that the process will settle with the point light source 1 10 at the desired distance from 
the SUT 140 (e.g., one focal length away). 

[000931 As noted above, the range-limiting aperture 160 clips out light that would otherwise 
produce an out-of-range condition for the wavefront senisor ISO. At the condition where the 
light is just beginning to clip, there is an additional effect due to diffraction around the edge 
of the aperture 160. This causes the image of the focal spot 188 to be non-symmetric at the 
detector array 1 84. 

[00094] This can be detected by analyzing an additional property of the focal spot 188 
incident on the detector array 184. The location of the focal spot 188 is usually determined 
by the application of a centroid or other spot locating algorithm. One algorithm is the first 
moment of the pixel values in the AOL This property (location) is compared against a 
reference location and used to determine the spot motion and hence the local wavefront slope. 
[00095] Another property that may readily be determined is the spot size. A change in the 
spot size, analogous to a change in position, may be used to determine if an aperture has 
affected the light that creates the spot. The spot size can be determined by additional analysis 
ofthe pixel values in the AOL One useful algorithm is the 2 moment algorithm. That is: 
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[00097] where each pixel is weighted by the square of the pixel position (relative to the 
center of the AOI), The signal Sy may be advantageously constructed from the actual pixel 
intensities lij by thresholding, deconvolution or other method to eliminate noise and crosstalk 
as appropriate. 

[00098] This algorithm produces a measure of the spot size in each direction (x and y). A 
measure of the skew or rotation can be obtained from the cross moment (x and y). 
[00099] This algorithm can be advantageously applied to identify a change in the spot size as 
compared to a reference. To this end, the spot size (in x and y, and xy as appropriate) is 
determined with a reference wavefront for each AOI. This is the same reference wave that is 
recorded when obtaining a reference set of spot locations for use in the wavefront slope 
calculation. The spot size is determined for each AOI using this (or other) methods for each 
frame of data that is analyzed. This is compared to the reference spot sizes. A difference 
(may be either larger or smaller) by more than a threshold amount may be the result of light 
cUpping on the edge of the aperture 160. Thus these AOIs have been identified, and may be 
removed computationally from fiirther calculations. This operates best in conjunction with 
the range linuted aperture 160, but may also be used by itself to identify and discard out-of- 
range AOIs. 

[000100] FIG. 4 shows an example AOI for a 4.6 mm focal length, 0.144 nam diameter 
lenslet array with 9.9 pm pixels. In this case the focal spot is well within range of the lenslet, 
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SO that it lies 100% enclosed in the boundaries of the AOL 

[000101] In contrast, FIG. 5 shows an AOI that has Ught from two different lenslets. This 
will clearly give a centroid location that is in error. This is the dangerous out-of-range 
condition which we seek to identify and avoid. It is dangerous because the centroid algorithm 
will produce a centroid location that is between the two focal spots which is completely 
wrong. The second moment in the y-direction for this case is 10 times the second moment of 
FIG. 4, making this easy to identify and eliminate with this method. 

[000102] Since the 2"^ moment algorithm works so well at identifying AOIs that have two 
(partial) focal spots, the range-limiting-aperture 160 may appear superfluous. This is not the 
case, however. If the local wavefront tilt is so great that the lenslet-to-AOI assignment is 
completely lost, such as when there is sufficient tilt for the AOI to contain the focal spot from 
the adjacent lenslet, the 2*^** moment would give no information. The centroid and 2^"^ 
moment look just like in FIG. 4, but the location of the focal spot is incorrect. The use of the 
range limiting aperture 160 prevents this from happening, since in this case the ligjht would be 
clipped by aperture 160 and the AOI would not contain any light that could lead to erroneous 
information. 

[000103] FIG. 6 illustrates such an intermediate case. In thus case, there is sufficient tilt that 
the focal spot is partially occluded by the range limiting aperture 160 (which was set to be 
slightly smaller than the AOI image in this instance). The focal spot has a distorted shape. In 
this case the 2"^^ moment is Vi the value in the y-direction as in the x-direction. This is readily 
detectable by comparison with the reference, and hence even these marginal AOIs can be 
identified and removed from the data that is used to compute the wavefront. 
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[000104] FIG. 7 shows the pattern of focal spots for a system such as that shown in FIG. 1 . 
In this case a very complex optic was used els the SUT 140. This particular optic had 
different optical power in 5 different rings. In this case the raijge-limiting aperture 160 served 
to eliminate the data where the different rings intercepted, The.local wave&ont slope at this 
point was multi-valued, so no good data could be obtained. The remaining data is all 
completely in-range, valid data and may be used to reconstruct the wavefront to achieve 
accurate results. 

[0001051 FIG. 8 is the same optic as tested in FIG. 7, but with a different location of the 
light source 110. In this case there was significant wavefront curvature. For constant 
curvature, the slope increases linearly with distance from the center of the optic. So at some 
radial point the wavefront slope exceeds the dynamic range of the wavefront sensor 180. In 
this case the range limiting aperture 160 has clipped all of this light so that only the central 
portion remains. While this does not coyer as large an area, it is more that sufficient to 
reconstruct the wavefront, determine the radius of curvature, and provide information 
feedback to set the location of the point light source 110. 

[000106] In the preceding discussion, it has been assumed that the point Ught source 1 10 is 
located on the optical axis of the SUT 140, and furthermore that the SUT 140 is aligned 
perpendicularly to the optical axis of the wavefront analysis system 100 as defined by the pair 
of lenses 150 and 170. In some applications, it is important that the point light source 1 10 is 
located fairly precisely on the optical axis of the SUT 140, and furthermore that the SUT 140 
is aligned fairly precisely perpendicularly to the optical axis of the wavefront analysis system 
100. 
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[000107] Accordingly, in one embodiment, the movable stage . 120 may be configured to not 
only move the point light source 110 aloi;ig the optical axis (z-aicis) of the wavefront analysis 
system 100, but also along the x-axis (side-to-side) and y-axis (up-and-down) directions with 
respect to the optical axis of the wavefront analysis system 100. Furthermore, the SUT holder 
130 may be configured to rotate and tilt the SUT 140 with respect to the optical axis of the 
wavefront analysis system 100. . . 

[000108] In that case, to adjust the point light source 110 until it is fairly precisely on-axis, 
the processor 190 moves the movable stage 120 between at least two different positions in the 
x-y directions, determines the location of the center of the light beam on the wavefront sensor 
180, and, from this data, calculates how far to move the point light source 1 10 to place it onto 
the optical axis of the wavefront analysis system 100. This procedure is typically performed 
with an aperture installed in the SUT holder 130 instead of an optical device or lens that is to 
be tested. The light beam passing through the aperture placed in the SUT holder 130 
produces a mmaiber of Ught spots (e.g., 20) on the wavefront sensor 1 80. By weighting each 
light spot and its associated location on the wavefront sensor^ 180, tiie center of the light beam 
can be determined, and the difference between the center of the light beam and the location 
and where the center should appear if the point light source 1 10 was on-axis (typically the 
center of the wavefront sensor 180) can be calculated. Alternatively, the beam center can be 
determined from the wavefront tilt detected by ttie wavefront sensor 1 80. This procedure can 
be performed in an automated manner with a closed loop feedback system including the 
processor 190 and the movable stage 120. 

[000109] Meanwhile, to fairly precisely align the SUT 140 to be perpendicular to the optical 
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axis of the wavefiront analysis system 100, the SXJT holder 130 may be rotated or tilted until 
the amount of astigmatism detected by the wavefiront sepsor 180 is minimized. The 
astigmatism maybe found, for example, frpm the coefficients of the Zemike polynomial for 
the detected wavefront. However, other outputs of the wavefiront sensor 180 maybe 
employed to determine when the SUT 140 is perpendicular to the optical axis of the 
wavefront analysis system 100. 

[000110] If the point light source 110 is located along the optical axis of the SXJT 140, you 
can only determine the on-axis aberrations. For most optical sj^tems, it is important to know 
both the on-axis and off-axis aberrations. This is especially important in camera lenses where 
good image quality is required at the comers of the image. By moving the point light source 
110 off-axis with the movable stage 120, the off-axis aberrations can be determined. Because 
of the imaging condition between the wavefiront sensor 180 and the SUT 140, the beam of 
light will remain on the wavefiront sensor 180 even when the point light source 1 1 0 is moved 
off axis. 

[000111] Furthermore, for some applications, it is also desired to know what effect the tip 
and tilt of the lens will have on the wavefiront aberrations. It is especially important to 
measure for the tip/tilt aberrations of intra-ocular lenses that ai'e put in the human eye after 
cataract surgery. Because of the large amoimt of variability of the structures that are foxmd in 
the eyes of different people, it is hard for the surgeons to control the tip and tilt of the intra- 
ocular Iots once it is implanted. So those kinds of lenses have to work well over a range of 
tip and tilt conditions. The tip/tilt wavefiront aberrations can be determined by adjusting the 
tip and tilt of the SUT 140 by means of the SUT holder 130 in- a determined pattem (e.g., 

29 



Docket No. WFS.OUC . .. 
-PATENT- 

under control of the processor 190) and making wavefront measurements at each tip and tilt. 
This SUT holder 130 may be adjusted manually, or under computer control with a motorized 
mount. 

[000112] When the point light source 1 10 is located one focal length ("f ) away from the 
SUT 140, the wavefront analysis system 100 measures the wavefront aberrations of the SUT 
140 for a coUimated light beam. However, many lenses, such ?is zoom lenses for example, 
are used with a variety of converging and diverging Ught beams, hi general, the aberrations 
of such a lens would change as the point light source 110 moves closer to or further from thp 
focal point of the lens (i.e., operate with convergmg or diverging light beams). Accordingly, 
it is sometimes desirable to characterize the aberrations of an SUT 140 with convergmg 
and/or diverging light beams. With the present system and method, the movable stage 120 
maybe moved along the optical axis (z-axis) of the wavefront analysis system 100 to pass 
converging or diverging light through the SUT 140, and the wavefront sensor 180 may make 
corresponding aberration measurements. 

[000113] When an intra-ocular.lens is used in an eye, it is immersed in fluid (primarily, 
water). The immersion changes the focal length of the lens to be about 2 to 3 times less than it 
is in air. In general the amount of .change depends on the index of refraction of the lens 
material, the thickness of the lens, and the front and back shapes of the lens. It is not 
possible to derive a general conversion formula from air to water unless those parameters are 
known. So, when those parameters are not known or have a great deal of uncertainty, it is 
very important to be able to measure the. lens when it is inmiersed in water. Accordingly, in 
one embodiment, the wavefix)nt analysis system 100 (and particularly the SUT holder 130) is 
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adapted to immerse the SUT 140 in water during the wavefront measurements. 
[000114] There is a new generation of multi-focal intra-ocular lenses that is being implanted 
in patients. These lenses are needed by many patients over the age of 40 that are presbyopic. 
Presbyopia is the condition where people can no longer adjust the focus of their eye between 
near and far. Presbyopia is a normal effect of aging. The traditional solution for these people 
is for them to wear glasses when they want to read, or for them to have bi-focal glasses. 
[000115] In the multi-focal intra-ocular lenses, the lens is designed so that different regions 
of the lenses have different focal lengths. The most common design is that the central 
circular region has one focal length, and the zone that is outside of that central region has a 
different focal length. In some designs, it might be that the top half of the lens has one focal 
length, and the bottom half has a different focal length. This is a very active area of research 
with a great many different types of designs being tested. 

[000116] The traditional technique for measuring these sorts of lenses is to physically place 
apertures in front of the lens and then do the measurenients. This is somewhat inconvenient 
as it requires making extra parts to measure the lenses, and for .tiiose parts to be properly 
aligned to the lenses. 

[000117] Accordingly, in one embodiment of the present invention, the processor 190 
operates with masks in the software that simulate the effect of physical apertures. Data in 
areas of the wavefront sensor 180 that is excluded by the mask(s) is not used when executing 
the algorithms described above for determining the focal.length. Thus, for example, the 
focal length of the central region of a multi- focal, intra-ocular lens can be first determined by 
masking out in software all of the wavefront data produced by the areas of the lens that are 
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outside the central region. Then, the focal length of the outer region of the multi-focal intra- 
ocular lens can be determined by masking out in software all of the wavefront data produced 
by the areas of the lens that are inside the central region. 

[000118] When the wavefront data from the wavefront sensor 1 80 is drawn on a display 
screen associated with the wavefront analysis system 100, the various regions can be outlined 
by a person using a mouse or pointer. But it is also possible for the processor 190 to execute 
a software algorithm to automatically determine where the masks should be, and for the 
analysis to be done separately over the different zones. Also, it possible to do optical analysis 
without the intermediate step of making internal masks, and instead making a zonal power 
map of the intra-ocular lens by using techniques that are similar to what are described in the 
Southwell paper, referenced above 

[000119] While preferred embodiments are disclosed herein, many variations are possible 
which remain within the concept and scope of the invention. For example, while the 
embodiment described above employs a movable platform to move the point light source with 
respect to a system or device under test, within a limited range (e.g., ±10 mm) it is possible 
to move the system or device under test, instead of the point light source, until the light beam 

emerging from the device or system is perfectly collimated,. Such variations would become. 

clear to one of ordinary skill in the art after inspection of the specification, drawings and 

claims herein. The mvention therefore is riot to be restricted ekcept within the spirit and 

scope of the appended claims. 
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